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Motivated by the recent experiment on kagome-lattice antiferromagnets, we study the zero-field 
ordering behavior of the antiferromagnetic classical Heisenberg model on a uniaxially distorted 
kagome lattice by Monte Carlo simulations. A first-order transition, which has no counterpart in 
the corresponding undistorted model, takes place at a very low temperature. Origin of the transition 
is ascribed to a cooperative proliferation of topological excitations inherent to the model. 



Renewed interest has recently been paid to novel order- 
ing properties of geometrically frustrated magnets [l|, Q ■ 
Two-dimensional (2D) kagome lattice, a corner-sharing 
network of regular triangles, is a typical example of such 
geometrically frustrated lattices. Many theoretical works 
have been performed to elucidate the ordering of the an- 
tiferromagnetic (AF) Heisenbergmodel on this lattice for 
both cases of quantum S = 1/2 [3|-l6| and classical S = oo 
[7 -11 1 spins. 



In the quantum S = 1/2 case, although there seems to 
be a consensus among researchers that the ground state 
of the model lacks in the conventional AF long-range or- 
der, there still exists a considerable diversity in the view 
about the nature of its ground state 0-@ • In the classi- 
cal S = 00 case, any combination of the local 120-degrees 
spin structure on a constituent triangle is allowed as its 
ground state, leading to the massive degeneracy. The 
effect of "order from disorder" is operative at low tem- 
peratures, which favors at the harmonic level the copla- 
nar states, or the spin nematic states, where all spins are 
contained on a common plane in the spin space [7|. Yet, 
the coplanar ground states themselves are heavily degen- 
erate, the degeneracy of 0(e N ) (N the number of spins). 
In order to determine which coplanar state is realized in 
the T — > limit, nonlinear excitations beyond the har- 
monic level need to be invoked. The model selects among 
them the so-called \/3 x \/Z state [1] . 

Experimentally, there have been few realizations of 
Heisenberg kagome AFs. Only recently, several possible 
realizations of S = 1/2 Heisenberg kagome AFs were re- 

II El, volborthite [IM! 



ported, e.g., herbertsmithite 

and vesigniete [n[. Volborthite Cu 3 V207(OH)2-2H 2 
and vesigniete BaCu3V20s(OH)2 are structurally dis- 
torted from the perfect kagome lattice in a uniaxial man- 
ner where the regular triangle is distorted to isosceles 
triangles. This structural distortion gives rise to two dis- 
tinct exchange couplings, J\ and J 2 , as demonstrated 
in the inset of Fig. 1(b), where the extent of the distor- 
tion may be represented by the parameter r = J 2 /Ji. 
Whether r is greater or smaller than unity in volborthite 
has not fully been determined experimentally, though the 
bond distance seems to favor r > 1 [HI]. An interesting 



recent experimental finding is that this compound ex- 
hibits a thermodynamic phase transition at T c ~ IK, 
with a slowly fluctuating ordered state EMI. 

Under such circumstances, it remains most interesting 
to clarify what type of ordering behavior is expected in 
the kagome Heisenberg AF under uniaxial distortion. In 
the present paper, we address this issue for the classical 
model in zero field. In particular, we wish to clarify what 
type of spin structure is chosen in the T — > limit, and 
whether a finite-temperature transition is ever possible. 

Only a few theoretical works have been made on the 
AF classical Heisenberg model on the distorted kagome 
lattice in the past. Y.B. Kim et al studied the model by 
an effective "chirality" Hamiltonian approach, and sug- 
gested that the so-called "chirality-stripe" state might be 
selected in the low-temperature limit [181 ] . Kancko et al 
studied the in-field properties of the model by Monte 
Carlo (MC) simulations, and observed a field-induced 
phase transition with a weak (zero)-field state being the 
\/3 x V3 state [HI]. 

In the present paper, we study the zero-field ordering 
properties of the model by means of extensive MC simula- 
tions, paying particular attention to the low-temperature 
regime which was not examined in previous works. Sur- 
prisingly, we find a clear first-order transition, which has 
no counterpart in the undistorted model, in the extremely 
low-temperature regime. 

The model we consider is the classical AF Heisenberg 
model on a uniaxially distorted kagome lattice with two 
distinct AF nearest-neighbor couplings J\ and J2 = rJi 
as shown in the inset of Fig. 1(b). The Hamiltonian is 
given by 



H 



Sj + J 2 £ Si 



(1) 



In the undistorted case r = 1 , the ground-state of a sin- 
gle triangle is the 120-degrees spin structure, whereas, in 
the distorted case r ^ 1, it is distorted such that a spin- 
canting angle at the Ji-bond becomes 9 = arccos(— ^-). 
The ground state of the distorted model is highly degen- 
erate in that any combination of the local ground state 
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on a triangle yields a ground state of the entire lattice 
fl8l [l9| . Order from disorder effect due to harmonic exci- 
tations then selects coplanar states as in the undistorted 
case. The coplanar states are still heavily degenerate, but 
the amount of degeneracy is much reduced from that in 
the undistorted case, from the extensive number of 0(e N ) 
to the subextensive number of 0(e L ) (L the linear dimen- 
sion) [18]. It should be noticed that the a/3 X V% state, 
which is the T — > state of the undistorted model, can- 
not be a ground state of the distorted model any more. 

In order to determine which state is favored in the 
T — > limit in the distorted model, we perform an ex- 
tensive MC simulation at low temperatures. Our MC 
simulation is a combination of the heat-bath and the 
over-relaxation methods. A unit MC sweep consists of 
single heat-bath sweep and subsequent 10 over-relaxation 
sweeps. Various lattice sizes are studied up to L = 192 
with total N = 3 x L 2 spins. We generate lxl 6 — 4x1 7 
MC steps per spin (MCS) at each temperature, both 
cooling and warming runs being made. For our check of 
equilibration, see the online Supplemental Material pol ]. 
We apply several types of boundary conditions (BC), in- 
cluding periodic, free, "fixed" and "vortex" boundaries 
(to be explained below). Most extensive calculations are 
performed for r = 1.1, while other values of r > 1 are 
also studied. 

In Fig. 1(a), we show the temperature dependence of 
the specific heat for the cooling runs for the case of 
r = 1.1. In the inset, the data are compared with those 
of the undistorted model r = 1. A sharp peak absent 
in the undistorted case appears, growing rapidly with L. 
The data turn out to be reversible for smaller sizes of 
L < 96, while a weak hysteresis is observed for larger 
sizes of L > 144. In Fig. 1(b), we show the temperature 
dependence of the energy for L = 192, where a clear 
hysteresis is observed in the transition region signaling 
the first-order nature of the transition. This is consistent 
with the sharp peak of the specific heat growing rapidly 
with L. The first-order transition is extremely weak in 
that the latent heat is only AE ~ 6 x 1CP 6 (in units 
of Ji). It is remarkable that such a small energy scale 
is generated in the model which has only O(l) energy 
scale. A closer look of the data reveales that the specific 
heat exhibits, in addition to the sharp main peak, addi- 
tional small peaks (or weak structures) at L-dcpendent 
temperatures. These small peaks are associated with a 
discontinuous shift of the g-value characterizing an in- 
commensurate spiral in finite systems under periodic BC, 
and are expected to vanish in the thermodynamic limit: 
See below. 

The inset of Fig. 2(a) exhibits the spin structure factor 
F in the (q x ,q y ) plane calculated at a temperature be- 
low the bulk trnasition temperature T c ~ 0.00065, where 
q is measured in units of 2ir/a (a the nearest- neighbor 
distance of the kagome lattice). The strongest inten- 
sity appears at a wavevector along the J2-bond direction 
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FIG. 1. (color online), (a) The temperature and size depen- 
dence of the specific heat per spin for the cooling run of the 
distorted model of r = 1.1 under periodic BC, either calcu- 
lated from the energy fluctuation or from the temperature 
difference of the energy. In the inset, the data for L = 60 
are compared with the corresponding data of the undistorted 
model r — 1 in a wider temperature range, (b) The tem- 
perature dependence of the energy of the distorted model of 
r = 1.1 for L = 192 both for the cooling and warming runs. 
The inset illustrates a distorted kagome lattice. 



slightly off the K point, corresponding to an incommen- 
surate spiral state. The observed peak has only a finite 
width reflecting a finite spin correlation length of the 2D 
Heisenberg model at finite T. Such an incommensurate 
spiral state is certainly a ground state of the distorted 
model, which reduces to the \/3 x v3 state in the undis- 
torted limit r — > 1. In the main panel of Fig. 2(a), we 
plot F(q x ,0) versus q x in the vicinity of the main peak 
for several temperatures across T c . Free BC are applied 
here in order to minimize the finite-size effect associated 
with the incommensurability As one approaches T c , the 
peak position gradually shifts from the K point to an 
incommensurate position. This shift of the peak posi- 
tion is observed already above T c for larger L. When 
one applies periodic BC, the ensuing discretization of q 
in units of 2tt/L often hampers the observation of such 
a shift for small L. This might explain the reason why 
the peak was located just at the K point in the simula- 
tion of Ref.[l9j]. For relatively large but finite L, such a 
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FIG. 2. (color online), (a) The inset exhibits the spin struc- 
ture factor in the wavevector (q x ,q y ) plane of the distorted 
model of r = 1.1 at a temperature below T c , T — 0.0005. 
Outer (inner) hexagon represents the first Brillouin zone of 
the triangular lattice of a lattice constant a (2a). Main 
panel plots F(q x ,0) versus q x in the vicinity of the K point 
q x /(2n) = 2/3 for several temperatures across T c . The lat- 
tice size is L = 96 with free BC. The vertical dotted lines 
represent the Appoint wavevector (right) and the incommen- 
surate wavevector at T = (left), (b) The temperature and 
size dependence of the nematic order parameter of the dis- 
torted model of r = 1.1 under periodic BC. In the inset, the 
L = 60 data are compared with the corresponding data of the 
undistorted model of r = 1. 



discretization effect of the periodic BC causes, with vary- 
ing the temperature, a sudden jump in the g-value. In 
any case, the low-temperature state of the model turns 
out to be an incommensurate spiral state, distinct from 
the chiral-spiral state inferred in Ref . fl8j . or from the 
\fi x V3 state suggested in Ref. [l9| . 

In Fig. 2(b), we show the temperature dependence of 
the nematic order parameter Q K representing the extent 
of the spin coplanarity 

(2) 

where N t is the number of upward triangles and the vec- 
tor chirality K a at each upward triangle a is given by 
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FIG. 3. (color online), (a) The temperature and size depen- 
dence of the number density of chiral domain walls of the 
distorted model of r = 1.1 under periodic BC. Inset exhibits 
a typical chiral-wall configuration. The symbols represent the 
upward and downward chirality vectors, (b) The temperature 
and the size dependence of the vortex energy AEtot(T) (the 
vortex free energy AFtot(T)) in the main panel (in the inset). 



k = t^(Si x 5 2 + 5 2 x S3 + S3 x Si), where three spins 

S\ ~ ^3 are arranged in a counter-clockwise manner on 
a triangle. 

In any 2D Heisenberg model, the nematic order pa- 
rameter should vanish at any finite T in the thermo- 
dynamic limit. Yet, it gives a measure of the nematic 
short-range order (SRO). As can be seen from Fig. 2(b), 
the nematic order or the spin coplanarity sets in almost 
simultaneously with the first-order transition for the case 
of r = 1.1. Additional step- wise structures of the data 
are due to the jump in the q- value associated with the in- 
commensurability effect in finite systems under periodic 
BC mentioned above. This onset of the nematic SRO 
occurs at a temperature significantly lower than that of 
the undistorted model: See the inset. 

Establishing the existence of a first-order transition, we 
further examine what is responsible for this transition. 
For this purpose, we study the behavior of two types of 
topological excitations inherent to the model. One is a 
chiral- domain wall, and the other is a Z2 vortex. 

The chiral domain wall is a wall-like excitation, which 
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might be obtained by connecting the reversal points of 
the vector chirality. It is a zero-mode arising from the 
coplanar nature of the low-T spin state. In the present 
distorted model, the chiral wall becomes a zero-mode 
only for a straight-line wall of system size L running per- 
pendicular to the J2-bond direction, in sharp contrast to 
the undistorted case where chiral walls of various sizes 
and shapes are possible including closed loops. 

Fig. 3(a) exhibits the temperature dependence of the 
total number density of chiral walls. Here we assign the 
vector chirality at each nearest- neighbor bond (ij) along 
the J2-bond direction as Hij = SiX Sj , and suppose that 
a lattice site i is included in the chiral wall when the con- 
dition /?i : i_|_e x • Ki_e X! i < is met. As the temperature is 
lowered across T c , the number of chiral walls diminishes 
rapidly. Typical snapshots of the chiral wall configura- 
tions are given in the online Supplemental Material plj ]. 

We also study the behavior of another topological exci- 



tation, Z2 vortex 22j , a point defect inherent to the frus- 



trated 2D Heisenberg model with the noncollinear spin 
order 

US-Si. 

For this purpose, we introduce two differ- 
ent BC: One is the fixed BC where we fix the boundary 
spins so that they match the expected ground-state spin 
configuration of an incommensurate spiral lying on the 
(S x ,S y ) plane. The other is the "vortex" BC where the 
boundary spins are fixed in the manner to accommodate 
a single vortex on top of the fixed-BC spin configura- 
tion, i.e., a spin rotation of 2tt around the S z -axis is 
applied to the boundary spins [26|. Then, we calculate 
the energy difference AE to t{T) between the total ener- 
gies of the fixed and the vortex BC, together with the 
total free-energy difference AF to t(T) evaluated by inte- 
grating AE tot (T) with respect to T 0. Fig.3(b) ex- 
hibits AE to t{T) calculated in this way. A sharp peak 
is observed just at T c , indicating that a strong anomaly 
occurs at T c in the vortex sector. The calculated free- 
energy difference AF to t(T) is shown in the inset, which 
indicates that the vortex tends to stiffen below T c . 

Thus, the transition turns out to accompany a signif- 
icant change in the state of topological excitations, i.e., 
accompany a proliferation (or a rapid increase in number) 
of both chiral walls and Z2 vortices. Indeed, as demon- 
strated in the online Supplemental Material |27[ , there is 
an effective attractive interaction acting between the chi- 
ral domain walls and the Z2 vortices, which might cause a 
cooperative thermal generation of these topological exci- 
tations possibly leading to the observed first-order tran- 
sition. The Z2-vortex transition often becomes of first 
order when it accompanies a simultaneous symmetry- 
breaking of other degrees of freedom 28, 29|. 

We have performed similar simulations for other val- 
ues of r > 1, and the resulting phase diagram is given 
in Fig. 4 in the distortion r vs. the temperature T plane. 
Together with the first-order transition line T c (r), the 
crossover line associated with the onset of the nematic 
order is also shown. Interestingly, this crossover line co- 
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FIG. 4. (color online). The phase diagram of the distorted 
kagome Heisenberg model in the temperature - distortion (r) 
plane. The undistorted model corresponds to r = 1. The 
crossover line is conventionally drawn at a point where Q K 
for L — 48 is equal to 0.01. The asymptotic behavior of 
the transition line near the undistorted limit is presently not 
clearly determined. 



incides with the first-order transition line for r > 1.1, 
while it is separated from the T c -line for r < 1.1, lying 
well above it. For smaller r < 1.04, we cannot identify 
the transition directly. This is either the transition it- 
self is absent for smaller r, or T c becomes so low that we 
simply cannot observe it in our simulations. In the range 
1.04 < r < 1.12, T c varies almost linearly with r. Lin- 
ear extrapolation yields a critical value r c ~ 1.02 below 
which wc no longer have a finitc-T transition. However, 
the difficulty to get reliable data at low temperatures in 
this small r region hampers us to reliable determine the 
fate of the transition line in the r — > 1 limit, i.e., whether 
i) it hits the T = axis at r = r c > 1, or ii) shows a non- 
linear behavior in the smaller r region reaching the T = 
axis only at r — 1, or iii) ends at a finite T exhibiting a 
critical end point at (r c > 1, T c > 0). 

Finally, we discuss the possible relation of the finite- 
T transition observed in our simulation to the zero-field 
transition of volborthite. The Curie temperature of the 
present model is estimated to be T cw ~ 0.9 for r = 1.1. 
If one matches it with the corresponding experimental 
value of volborthite ~ 115K [30j |. T c is estimated to be 
0.1K, which is an order of magnitude smaller than the 
experimental value, IK. Several possible causes for this 
discrepancy are conceivable: i) The experimental value 
of the distortion r is largely unknown which could ex- 
plain at least a part of the discrepancy, ii) Volborthite is 
a S = 1/2 system. Quantum effect neglected here might 
give significant correction, iii) Real material is 3D. Weak 
interplane coupling neglected here might push up the T c - 
value significantly, iv) The Dzyaloshinskii-Moriya inter- 
action neglected here might change the ordering behavior 
significantly. Further studies are then required to clarify 
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the true relation between the transition found here and 
the experimental one. 

In summary, we studied the ordering of the classical 
Heisenberg AF on a distorted kagome lattice by MC sim- 
ulations. We have found that the model exhibits a ther- 
modynamic first-order transition at an extremely low but 
finite temperature. Cooperative generation of the chiral 
domain walls and the Zi vortices might be responsible 
for this first-order transition. 
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Supplemental Material 



In this Supplemental Material, we explain some of the 
details about (I) the equilibration check of our Monte 
Carlo simulations, (II) the states of the chiral domain 
walls, (III) the properties of the Z 2 vortices, and (IV) 
the manner how the Z2 vortices and the chiral domain 
walls interact with each other in our model. 



EQUILIBRATION CHECK OF MONTE CARLO 
SIMULATIONS 



In this section, we explain how we check equilibration 
in our Monte Carlo (MC) simulations. Since our simula- 
tions are performed at very low temperatures, an appro- 
priate check of equilibration is crucially important. The 
standard check might be to monitor the stability of vari- 
ous observables as a function of MC time for successively 
longer observation times, which we do perform. In addi- 
tion to this standard test, we also perform the following 
procedure for an additional check of equilibration. 
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FIG. 5. The temperature dependence of the specific heat cal- 
culated either by the heat-bath (HB) method combined with 
the over-relaxation method, or by the Wang-Landau (WL) 
method. The distortion parameter is r = 1.1. 

As mentioned in the main text, we take our data based 
on the combination of the heat-bath (HB) method and 
the over-relaxation method. (The temperature-exchange 
technique often useful in, e.g., spin-glass problems, is not 
very efficient here at such low temperatures.) To be sure 
that our data are fully thcrmalized, we also take the data 
of the same physical quantities based on the so-called 
Wang-Landau method jl[ and check whether the data 
taken by the two different methods agree within the er- 
rors. Note that the Wang-Landau method directly com- 
putes the density of states, and the temperature comes 
into physical quantities in a way very different from the 
heat-bath (+over-relaxation) method. Thus, the agree- 
ment of the datasets taken by these two different methods 



should yield a stringent test of equilibration. In Fig. 5, we 
show the temperature dependence of the specific heat for 
the case of r = 1.1 calculated either by the heat-bath (+ 
over-relaxation) method or by the Wang-Landau method 
for the sizes L = 12 and 24. As can be seen from the fig- 
ure, the two kinds of datasets agree completely, demon- 
strating that our data are a fully thermalized even at the 
lowest temperature studied. For larger sizes, the Wang- 
Landau method is very much time-consuming, and we 
are unable to obtain the data for larger sizes. Yet, even 
for smaller sizes of L — 12 and 24 shown here, the exis- 
tence of the growing specific-heat peak is already clear, 
demonstrating the occurrence of an equilibrium phase 
transition at this very low temperature. 



CHIRAL DOMAIN WALLS 

In this section, we show how the chiral domain walls 
are thermally generated in our model both above and 
below T c . In Fig. 6, we give typical snapshots of the chi- 
ral domain-wall configurations below T c , T = 5 x 10 -4 
[left figure] and above T c , T — 7 x 10 -4 [upper figure] 
for a system of the distortion parameter r = 1.1 where 
T c ~ 6.5 x 10~ 4 . Free boundary conditions (BC) are 
applied here, because in general incommensurate cases 
periodic BC applied to finite lattices force the system 
to introduce a pair of chiral domain walls which are ex- 
tended in the perpendicular direction. By contrast, no 
such constraint exists in free BC which we adopt here. 
Chiral domain walls, defined at a lattice point i connect- 
ing two bonds (i — e^,i) and (i,i + e^) in line in the 
manner as explained in the main text, are depicted in 
the figure as blue segments perpendicular the bond and 
passing the site i. 

Chiral domain walls may be classified into two types: 
An extended wall spanning the entire lattice from one 
end to the other perpendicular to the J2-bond direction 
and a closed-loop wall. Note that a closed-loop wall costs 
a small but nonzero energy, while a perfectly straight-line 
wall perpendicular to the J2-bond direction costs no en- 
ergy (zero mode) owing to the large degeneracy of the 
ground state of the J\ — J2 kagome Heisenberg model. 
As can be seen from the figure, both types of chiral do- 
main walls are present below and above T c , whereas their 
number density increases rapidly as the temperature is 
increased across T c . As can be seen from the upper fig- 
ure, an extended wall running along the perpendicular 
direction seems to persist in small number even below T c 
(note that free BC are applied here so that the observed 
chiral domain wall is not the one forced by periodic BC). 
It may be no wonder that extended walls remain in small 
number even at low temperatures, since they are zero- 
energy modes so long as they run straight along the per- 
pendicular direction. 
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FIG. 6. Typical snapshots of the chiral domain-wall con- 
figurations both below T c , T = 5 x 10~ 4 [upper figure], and 
above T c ,T = 7xlCP 4 [lower figure]. The distortion param- 
eter is r = 1.1. Free BC are applied for a system of linear size 
L = 96. 

Z 2 VORTICES 

In this section, we explain some of the basic prop- 
erties of the Z2 vortex, together with the way how we 
identify the spatial position of the Z2 vortex (more pre- 
cisely, the position of its core) in our simulations. The 
Z-2 vortex is a topologically stable point defect inherent 
to two-dimensional frustrated Heisenberg systems with 
the noncollinear spin order. It possesses a parity-like Zi 
(±1) topological quantum number depending on whether 
there is a vortex or no vortex. As such, a winding num- 
ber, familiar as a good topological quantum number char- 
acterizing the standard vortex, is not a good topological 
quantum number here. 

The properties of the Zi vortices were discussed in 
some detail in Ref. [2] in the context of the classi- 
cal Heisenberg antiferromagnet (AF) on the triangular- 
lattice with the nearest-neighbor coupling. The Z% vor- 
tex might be regarded as a vortex formed by the chirality 
vector. It was suggested in Refs.0, HJ that the triangu- 
lar Heisenberg AF might exhibit a thermodynamic phase 
transition at a finite temperature driven by the binding- 



unbinding of the Z2 vortices, keeping the standard spin 
correlation length to be finite. In the previous studies, 
as an order parameter characterizing the ^-vortex tran- 
sition, either a Wilson-loop (a vorticity function) Q or a 
vorticity modulus 0, Q has been proposed. The vortic- 
ity modulus v might be defined as the total free-energy 
cost against a vortex formation AF to t divided by hiL, 
v = AF tot /lnL. In the triangular- lattice Heisenberg 
AF in its thermodynamic limit, v should become zero in 
the high-temperature phase where there exist free vor- 
tices, while v should take a nonzero value in the low- 
temperature phase where all vortices are paired with no 
free vortices. In case of our kagome model, the existence 
of the chiral domain-wall excitations would modify the 
nature of the vortex transition as wc shall see in the next 
section. 

In our simulations on a (distorted) kagome Heisenberg 
AF, we monitor the spatial position of the Z2 vortex fol- 
lowing the procedures of Ref.0]. The triangular lattice 
consists of three interpenetrating triangular sublatticcs, 
and one might take an elementary plaquette a minimum 
upward triangle consisting of the three sublattice sites. 
At low temperatures, the local spin structure at each 
plaquette is approximately a 120-degrees structure. In 
Ref. [2j, an 5*0(3) local 'frame' was assigned to each such 
elementary plaquette, and a Z2 vortex is identified for 
every minimum triangular loop consisting of three such 
plaquettes. We modify such a procedure of the local vor- 
tex identification to our kagome system. 

First, an elementary plaquette is taken again as a min- 
imum upward triangle consisting of three spins, which 
coincides with a crystallographic unit cell of the kagome 
lattice. The lattice formed by these elementary plaque- 
ttes also forms the triangular super-lattice consisting of 
the three triangular sublatticcs A, B and C. The spin 
structure here is an incommensurate spiral state close to 
the y/3 x \/3 state, as demonstrated in the main text. 
In the perfect -\/3 x state, the 120-degrees structure 
on an elementary plaquette would exhibit a uniform ar- 
rangement on each sublattice A, B and C. In view of this, 
we assign our local 'frame' to an elementary plaquette on 
one of these three sublattices, say, on A, and identify an 
elementary Z2 vortex for every minimum triangular loop 
consisting of three such plaquettes of the sublattice A. 
While the true spin structure generally deviates from the 
y/3 x s/3 structure, this deviation hardly affects the defi- 
nition of the Z2 vorticity since only the topological char- 
acter matters in the definition. In computing the spatial 
distribution of the Z2 vortices in the following section, 
we use this definition of the local Z2 vorticity. 

Z 2 VORTICES AND CHIRAL DOMAIN WALLS 

In this section, we explain the manner how the Z2 vor- 
tices and the chiral domain walls interact with each other 
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FIG. 7. The temperature dependence of the number density 
(the appearance probability) of Z2 vortices for the two kinds 
of sites, the site attached to the chiral domain wall and the site 
away from the chiral domain wall. The distortion parameter 
is r — 1.1, periodic BC being applied. The inset represents 
the temperature dependence of the number density of the 
Z2 vortices for the two kinds of chiral domain-wall sites, the 
site attached to an extended, system-spanning wall, and the 
site attached to a closed-loop wall. The arrow indicates the 
transition point. 



both above and below T c . 

In Fig. 7, we show the temperature dependence of the 
number density (or the appearance probability) of Z2 
vortices, d v , for the two kinds of sites, the site attached 
to the chiral domain wall and the site away from the chi- 
ral domain wall. The distortion parameter is set r = 1.1. 
As can immediately be seen from the figure, the sites at- 
tached to the chiral domain wall allow higher density of 
Z2 vortices. Namely, the Z2 vortices are generated pri- 
marily at or near the chiral domain walls. It means that 
there exists an effective attractive interaction between the 
Z2 vortices and the chiral domain walls. The inset repre- 
sents the temperature dependence of the number density 
of Z2 vortices for the two kinds of chiral domain-wall 
sites, the site attached to an extended, system-spanning 
domain wall and the site attached to a closed-loop wall. 
One can see from the inset that at lower temperatures an 
extended chiral domain wall accommodates a Z2 vortex 
preferably. 

We note that the existence of an effective attractive in- 
teraction acting between the chiral domain walls and the 
Z2 vortices might provide a possible explanation of the 
observed first-order transition. As the temperatures is in- 
creased, both the chiral domain walls and the Z2 vortices 
tend to be thermally excited. If there is a strong attrac- 
tive interaction between these two topological excitations 
as we have demonstrated above, one would expect a pos- 
itive feed-back effect operative in their thermal genera- 
tion: i.e., the generation of one induces the generation of 
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FIG. 8. Snapshots of the ^2-vortex configurations both be- 
low T c , T = 5 x 1CT 4 [top and middle fi gures], and above T c , 
T = 7 x 10~ 4 [bottom figure], depicted by black dots. Chiral 
domain wall are also shown by blue lines. The distortion pa- 
rameter is r = 1.1, periodic BC being applied. All figures are 
parts of the L — 96 lattice. The top figure is a part of the lat- 
tice containing an extended chiral domain wall, whereas the 
middle figure is another part of the same lattice not containing 
an extended chiral domain wall, i.e., containing closed loops 
only. Below T c , while an extended wall allows isolated, ap- 
parently free Z2 vortices, Z2 vortices appear always in bound 
pairs at closed-loop wall sites. 



the other. Thus, the observed first-order transition might 
be a consequence of such cooperative thermal generation 
of the chiral domain walls and the Z2 vortices. 

In Fig. 8, we show snapshots of the Z2-vortex config- 
urations both below T c , T = 5 x 10 -4 [top and middle 
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figures], and above T c , T = 7 x [bottom figure], 

together with the corresponding chiral domain-wall con- 
figurations. The distortion parameter is set r = 1.1. All 
figures are parts of the L = 96 lattice. The top and the 
middle figures represent different parts of the same lattice 
measured at exactly the same MC time. The top figure is 
a part of the lattice containing an extended chiral domain 
wall, whereas the middle figure is a part not containing 
an extended chiral domain wall, i.e., a part containing 
closed-loop walls only. As can be seen from the top fig- 
ure, an extended wall allows isolated, apparently free Z2 
vortices even below T c in addition to ^2-vortex pairs, 
whereas closed-loop walls bear only i^-vortex pairs be- 
low T c . Above T c , as can seen from the bottom figure, 
free vortices are generated not only at extended- wall sites 
but also at closed-loop wall sites, whereas the number of 



Z2 vortices becomes large and the distinction between 
free vortices and vortex-pairs become somewhat obscure 
here. 
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